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ABSTRACT 

Solutions  are  given  for  the  partial  differential  equation 
-Jjr  u(t,x,y)  +  ■—  f  (u  (t  ,x, y ) )  +  g(u(t,x,y))  =  0 

with  initial  data  constant  in  each  quadrant  of  the  (x,y)  plane.  This  problem 
generalizes  the  Riemann  Problem  for  equations  in  one  space  dimension.  Although 
existence  and  uniqueness  of  solution  are  known,  little  is  known  concerning  the 
qualitative  behavior  of  solutions. 

When  f  and  g  are  convex  and  f  =  g,  then  our  solutions  satisfy  the 
uniqueness,  or  entropy  condition  given  by  Kruzkov  and  Vol'pert.  Under  certain 
extra  conditions  on  f  and  g,  our  solutions  satisfy  the  entropy  condition 
if  f  and  g  are  convex  and  sufficiently  close.  A  counterexample  is  given 
to  show  the  necessity  of  these  extra  conditions  on  f  and  g.  The  correct 
entropy  solution  for  this  counter-example  exhibits  new  and  interesting  phenomena 
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SIGNIFICANCE  AND  EXPLANATION 


In  making  numerical  computations  of  multi-dimensional  gas-dynamic  flows 
involving  shock  waves,  a  major  difficulty  arises  when  two  shock  waves  cross 
each  other.  We  have  made  a  modest  beginning  on  problems  of  this  type,  in  the 
case  of  a  small  class  of  single  conservation  laws,  by  giving  exact,  explicit 
solutions  of  two  dimensional  Riemann  Problems.  Our  initial  data  consist  of 
step  functions  constant  in  the  quadrants  of  the  (x,y)  plane.  However,  our 
results  apply  to  any  two  intersecting  lines  of  discontinuity. 

Our  results  may  also  be  potentially  applicable  to  the  study  of  oil  recovery 
problems,  with  regard  to  "fingering"  phenomena. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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THE  RIEMANN  PROBLEM  IN  TWO  SPACE  DIMENSIONS 
FOR  A  SINGLE  CONSERVATION  LAW 

* 

David  H.  Wagner 

1.  Introduction. 

Let  f  and  g  be  given  real  functions  satisfying  f"  >  0  and  g"  >  0. 
Consider  the  initial  value  problem 

■|j:  u(t,x,y)  +  f (u(t,x,y))  +  ^  g(u(t,x,y))  »  0,  (1.1) 

u(0,x,y)  =  u^  for  x  >  0,  y  >  0  ,  (1.2) 

*  Uj  for  x  <  0,  y  >  0  , 

*  for  x  <  0,  y  <  0  , 

=  for  x  >  0,  y  <  0  . 

This  is  a  Riemann  Problem  in  two  space  variables.  It  generalizes  the 
Riemann  Problem  in  one  space  variable,  the  study  of  which  has  been  a 
key  to  the  understanding  of  solutions  to  systems  of  nonlinear  hyperbolic 
conservation  laws  in  one  space  variable  [2). 

Global  existence  of  weak  solutions  to  (1.1)  with  more  general 
initial  data  than  (1.2)  was  first  proved  by  Conway  and  Smoller  [1]. 

Later  Vol'pert  [9],  and  Kruzkov  [6]  proved  existence  and  uniqueness  of 
weak  solutions  satisfying  an  entropy  condition,  in  the  class  of  bounded 
measureable  functions. 

No  similar  advances  have  been  made  concerning  systems  of  non¬ 
linear  hyperbolic  conservation  laws  in  two  or  more  space  variables;  it 
may  be  that  study  of  the  Riemann  Problem  for  these  systems  will  yield 
a  breakthrough.  In  this  paper  we  begin  an  attack  on  this  problem  by 
finding  explicit  entropy  solutions  to  (1.1),  (1.2)  for  a  la^ge  class 
of  pairs  (f,g),  in  the  case  of  a  scalar  conservation  law. 

We  should  mention  that  Guckenheimer  [5],  and  Val’ka  [8],  have 
*  - 
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studied  examples  of  (1.1)  with  piecewise  constant  initial  data,  in  con¬ 
figurations  different  from  (1.2). 


+  2 

Definition  1.1.  A  bounded,  measurable  function  u:  1R  *  TR  ■*  1R 
is  said  to  be  a  weak  solution  to  the  initial  value  problem  consisting  of 
(1.1)  with  initial  data  u(0,x,y)  =  Uq(x,y)  if 

j^+  ^  u  ^  +  f(u)  ^  +  g(u)  ^  $  dx  dy  dt  -  0,  (1.3) 

,  OO  -f  2  V  1 

for  every  test  function  $  e  Cq(1R  x  1R  ),  and  if  u(t,*,*)  -►  in  L^oc 

as  t  +  0. 


Definition  1.2.  (Vol'pert  [8],  Kruzkov  [6]).  A  weak  solution 
u,  to  (1.1),  is  said  to  satisfy  the  entropy  condition  if  for  any 
real  constant  k  and  any  e  C*(1R+  x  1R2)  such  that  §  >  0, 

|  +  |  2  sign(u-k)[(u-k)  ^  $  +  (f(u)  -  f(k))  (j)  (1.4) 

1R  1R 

+  (g(u)  -  g(k))  ~  <)]  dx  dy  dt  S  0  . 

We  shall  construct  weak  solutions  to  (1.1),  (1.2)  which  are  valid 
for  any  choice  of  the  constants  u^,  u^,  u^,  and  u^.  It  is  conceivable 
that  this  construction  may  fail  to  produce  a  well  defined  function  if 
f  and  g  are  not  sufficiently  close  to  each  other.  The  degree  of  close¬ 
ness  that  is  sufficient  may  depend  on  the  choice  of  f  and  g;  therefore, 
in  order  to  be  rigorous  we  will  state  and  prove  our  theorems  in  terms 
of  the  distance  from  f  and  g  to  a  given  reference  function  h.  Thus  we 
consider  the  pair  (f,g)  as  a  perturbation  from  the  pair  (h,h).  Of  course, 
our  theorems  will  cover  the  case  where  f  is  held  fixed  and  g  is  perturbed 
away  from  f. 
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The  form  of  our  solution  to  (1.1),  (1.2)  varies  with  different 
orderings  of  the  constants  u^.  Thus  there  are  twenty-four  cases  to  be 
considered.  Fortunately,  these  twenty-four  cases  can  be  reduced,  via 
geometrical  reflections  and  inversions,  to  eight.  We  will  however, 
consider  each  of  the  twenty-four  cases  and  Identify  the  reductions. 

We  shall  show  that  our  construction  produces  a  well  defined  function 
which  satisfies  the  entropy  condition,  if  f  =  g,  and  f"  >  0;  see  Theorem 
1.  We  will  also  show,  in  Theorem  2,  that  under  certain  ordering  condi¬ 
tions  on  u^,  u^,  u^,  and  u^,  we  have  that  for  every  function  h  such  that 

h"  >  0,  there  exists  e  >  0  such  that  [|f-hj|  -  <  £  and  ||g-h||  ?  <  e 

C  C 

imply  that  our  construction  produces  a  well  defined  function  which  satisfies 

the  entropy  condition.  In  Theorem  3  we  shall  show  that  if  the  ordering 

conditions  of  Theorem  2  are  not  satisfied,  then,  provided  f"  =  g"  and 

f'"  =  g'"  at  certain  points  w  to  be  specified  later,  for  every  h  such 

that  h"  >  0,  and  h"  =  f"  and  h"'  =  f"'  at  all  of  the  points  w,  there  exists 

e  >  0  such  that  ||  f  —  h i J  ,  <  e  and  ||g-h||  ,  <  e  imply  that  our  con- 

C  C 

struction  produces  a  well  defined  function  which  satisfies  the  entropy 
condition. 

Finally,  an  example  will  be  given  where  f'"  ^  g'M  at  the  point  w  of 
Theorem  3,  and  where  our  constructed  solution,  although  it  is  a  function, 
does  not  satisfy  the  entropy  condition.  In  this  example  f  and  g  may  be 
arbitrarily  close,  and  the  initial  data  may  be  arbitrarily  small.  We 
will  also  give  the  correct  entropy  solution  for  this  example. 

One  system  of  equations  containing  a  scalar  conservation  law  is  that 
describing  two  phase,  two  dimensional  immiscible  flow  in  porous  media. 
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where  gravity,  capillary  pressure,  molecular  diffusion,  compre:  ibllity, 
as  well  as  spatial  variations  in  porosity,  depth,  and  viscosity  have  been 
neglected  (see  (3,4,71): 


(vlf(S))  +  ^  (v2f(S))  =  0  (1.3) 

v  =  (v1»v2>  =  -k(S)7p  (1.6) 

-*  -»■ 

V  •  v  =  source  terms  .  (1.7) 


2 

In  r 3,4  I,  f(S)  =  s  /k(S)  was  used.  Although  one  may  imagine  that  our 
solutions  are  thus  special  entropy  solutions  of  this  system  with  Vp  con¬ 
stant,  (1.7)  prevents  this.  However  some  of  our  solutions,  Cases  9  and 
19, exhibit  shock  waves  meeting  at  an  accute  angle,  or  in  a  cusp,  similar 
to  the  fingering  behavior  which  is  of  interest  in  oil  recovery  problems. 


and  for  which  (1.5  -  1.7)  is  a  model. 


2.  Construction  of  the  Solutions. 

We  shall  see  that  our  constructed  solutions  are  piecewise  smooth, 
having  discontinuity  sets  consisting  almost  everywhere  with  respect  to 
two  dimensional  Hausdorff  measure,  of  smooth  surfaces.  In  this  context, 
Definition  1.2  implies  two  conditions,  given  below,  on  the  discontinuities 
of  a  solution.  These  can  easily  be  derived  via  localization  and  integra¬ 
tion  by  parts,  and  appropriate  choices  of  the  constant  k. 

Condition  2.1.  (The  Rankine-Hugoniot  condition).  At  any  point 
p  on  a  surface  of  discontinuity  S  of  the  solution  u,  if 

(a)  n  is  a  unit  normal  vector  to  S  at  p, 

(b)  u+  =  lim  £  ->  0+  u (p  +  en)  , 

(c)  u  =  lim  e  -*■  0+  u(p  -  en)  , 
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then 

n  •  (u+  -  u  ,  f(u+)  -  f (u  ),  g(u+)  -  g(u  ))  =  0  .  (2.1) 

-y  -4"  — 

Condition  2.2.  (The  entropy  condition) .  Orient  n  so  that  u  >  a  . 

If  k  is  any  constant  such  that  u  'S  k  5  u"*",  then 

n  •  (k  -  u+,  f(k)  -  f(u+),  g(k)  -  g(u+))  >  0  .  (2.2) 

Using  (2.1)  one  may  check  that  (2.2)  is  equivalent  to 

n  •  (k  -  u  ,  f (k)  -  f(u~),  g(k)  -  g(u  ))  >  0  .  (2.3) 

One  may  further  check  that  if  a  function  u  is  a  piecewise  classical  solution, 
except  for  smooth  surfaces  of  discontinuity  where  Conditions  2.1  and  2.2 
hold,  then  u  is  a  weak  solution  satisfying  the  entropy  condition. 

Let  us  consider  one  dimensional  shock  waves  and  rarefaction  waves, 
as  they  arise  in  the  two  dimensional  Riemann  Problem. 

(a)  The  one  dimensional  shock  wave.  If  the  initial  data  is 

u(0,x,y)  =  Uj  if  x  <  0  ,  (2.4) 

=  U£  if  x  >  0  , 

and  Uj  >  \x  then  the  problem  really  has  only  one  space  dimension: 

ut  +  f  ^x  =  0  ’  (2.5) 

u(0,x)  =  u^  if  x  <  0  , 

=  u^  if  x  >  0. 

In  this  case  the  solution  is  well  known: 
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This  solution  has  a  discontinuity,  called  a  "shock  wave,"  along  the  plane 

(2.7) 


f(uj)  -  f(u2) 
x  = - t  . 


U1  '  U2 


We  will  refer  to  this  shock  wave  as  "SXtLupUjU,"  for  "shock  in  the  x 
direction  connecting  u^  to  u2."  The  shock  wave  in  the  y  direction  obtained 
by  interchanging  x  with  y  and  f  with  g  above  we  will  call  "SYEu^ ,u2 J . " 

(b)  The  one  dimensional  rarefaction  wave.  If  the  initial  data  is: 


u(0,x,y)  =  Uj,  if  x  <  0  ,  (2.8) 

“  u2>  if  x  >  0  , 

and  Uj  <  u^,  then  the  problem  again  reduces  to  (2.5),  and  its  solution  is 
well  known: 


u(t,x,y)  =  u(t,x)  =  Uj  if  x  <  f’^H,  (2.9) 

■  Uj  if  x  >  f'  (u2)t  , 

=  8  if  X  =  g'(s)t  ,  U^  ^  s  S  U2  . 

The  part  of  this  solution  between  x  =  f'(u^)t  and  x  =  f'(u2)  is  called  a 
"rarefaction  wave."  We  will  refer  to  this  particular  rarefaction  wave  as 
RX[Uj,u2],  for  "rarefaction  in  the  x-direction,  connecting  u^  to  u2." 

Note  that  the  solution  described  in  (a)  is  a  weak  solution  to  (b) , 
since  it  satisfies  Condition  (2.1).  However  it  does  not  satisfy  the 
entropy  condition  since  u^  <  u2>  and  f  is  assumed  to  be  convex. 


The  interaction  of  RX  with  SY.  Let  the  initial  data  be  as  in 
(1.2),  with  u^  ■  u2  ■  u^  *  w,  and  u^  *  v,  and  v  >  w.  Then  for  bounded 
t,  and  sufficiently  large  x,  the  solution  looks  locally  like  SY[v,w], 


due  to  the  principle  of  finite  domain  of  dependence,  which  was  shown  to 
hold  in  this  context  by  Vol'pert  [9],  and  Kruzkov  [6],  For  y  sufficiently 
negative,  the  solution  looks  locally  like  RX[w,vJ.  Since  a  solution  is 
invariant  under  dilations  (t,x,y)  -*  (ct.cx.cy)  for  c  >  0  whenever  the 
initial  data  u(0,x,y)  is  invariant  under  dilations  (x,y)  •+■  (cx,cy),  we 
may  describe  a  solution  completely  by  describing  it  along  the  plane  t  =  1. 
The  solution  is  constant  on  rays  through  the  (t,x,y)  origin. 

Thus  our  current  knowledge  of  the  solution  to  this  problem  may  be 
described  by  Figure  1.  In  Figure  1,  the  horizontal  line 


y 


RX[w,v] 

Figure  1.  Interaction  of  RX  with  SY 


labelled  "SY[v,wJ"  indicates  the  plane  of  that  shock  wave,  and  the  vertical 
lines  labelled  RX[w,vl  indicate  planes  along  which  the  solution  u  is  con¬ 
stant,  thereby  depicting  a  rarefaction  wave.  The  space  labelled  "?"  is 
filled  in  as  follows.  The  RX  region  meets  the  region  were  u  =  w  along  a 
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smooth  surface  of  discontinuity  S,  having  equations  x  =  f'(s)t,  v  =  v(s)t, 
parameterized  by  s,  t  for  w  <  s  t  v.  The  unknown  function  y  is  determined 
by  the  jump  conditions,  is  follows. 

First  we  describe  the  normal  vector  n  =  (n  ,n  ,n  ),  to  the  shock 

t  x  t 

surface  s  in  terms  of  f  and  y.  To  do  this,  we  need  two  tangent  vectors 
to  the  surface  x  =  f'(s)t,  y  =  y(s)t. 

Holding  s  fixed,  we  have  dx  =  f'(s)dt,  dy  =  y(s)dt.  Holding  t 
fixed,  we  have  dx  =  f"(s)tds  ,  dy  =  y'(s)tds  .  Thus  we  have  two  tangent 
vectors,  ^  =  (1,  f'(s),  y(s)),  v2  =  (0,  f"(s),  y'(s)).  Inen 

n  =  vj  x  v2  -  (-f"(s)y(s)  +  f'(s)y’(s),  -y’(s),  f"(s)).  (2.10) 

Keeping  in  mind  that  in  RX[w,v|,  u  =  s  on  the  plane  x  =  f'(s)t  for 
y  <  y(s)t,  the  Ranklne-Hugoniot  condition  gives  us: 

(W  -  s)(f'(s)y'(s)  -  f"(s)y(s>)  +  (f (w)  -  f (s)) (-y ' (s))  (2.11) 

+  f"(s)(g(w)  -  g (s) )  =  0  . 


This  equation  is  a  first  order,  linear,  scalar  differential  equation 
for  the  unknown  function  y: 


y’(s) 


f../_xg(w)  -  g(s)  -  y(s)(w  -  s)  . 
V  'f(w)  -  f  (s)  -  f ’ (s)  (w  -  s) 


(2.12) 


Note  that  the  denominator  in  the  right  hand  side  of  (2.12)  is  always 
positive  for  w  ^  s,  since  f"  >  0  . 

The  shock  surface  S  should  be  a  smooth  continuation  of  the  planar 
shock  wave  SY[v,w]  through  the  rarefaction  wave  RX[w,vl.  Therefore  these 
two  surfaces  should  meet  along  a  common  line  at  the  right  edge  of  RX[w,vi. 
This  yields  the  following  initial  condition  for  y: 
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Y(v)  = 

'  V  -  w 


(2.13) 


Since  (2.12)  Is  linear,  it  is  thus  explicitly  solvable.  The  following 
formula  for  y  can  be  obtained  from  the  usual  one  by  an  integration  by  parts: 


y  (s)  =  B(w)  (2.H) 

s  -  w 

f"(z)(w  -  2) _ 

f  (w)  -  f  (z)  -  f '  (z)  (w  -  z) 

.  g(w)  -  g(r)  -  g'  (r)  (w  -  r)  dj_ 

(w  -  r)2 

Note  that  for  w  <  s  <  v,  we  have  Y(s)  >  ^ ^ s-- — —  B(v^  ;  hence, 

s  -  w 

using  (2.12),  we  have  y'(s)  >  0.  Therefore,  as  s  decreases  to  w,  y(s) 
approaches  some  finite  limit,  which  is  greater  than  or  equal  to  g' (w) . 

In  fact  this  limit  is  g'(w),  as  we  shall  see  in  §3. 

We  shall  refer  to  the  shock  surface  S  as  "r[v,w]."  See  Figure  2. 

This  shock  surface  will  occur  in  solutions  to  the  Riemann  Problem,  and 
usually  in  truncated  form,  denoted  here  r[v,w:p],  namely  the  portion  of 
rCv.wl  corresponding  to  values  of  s  greater  than  some  number 
p  and  less  than  v.  Reflecting  T[v,w]  across  the  line  x  =  y  by  inter¬ 
changing  x  with  y  and  f  with  g,  we  obtain  a  similar  shock  wave,  denoted 
'TRrv.w]." 


Note  that,  given  any  solution  u  to  (1.1),  with  initial  data 
uQ(x,y),  the  function  u  given  by  u(t,x,y)  =  -u(t,-x,-v)  is  a  solution  to 


ut  +  f(-u)x  +  g(-u)y  =  0 

with  initial  data  u(0,x,y)  =  -U()(-x,-y).  Thus,  if  u(t,x,y)  = 

I’(t,x,y, f .g.Uj.u^.u^.u^)  is  a  formula  giving  the  solution  to  the  Riemann 


Figure  3.  The  Shock  Surface  n[v,w] 


We  are  now  ready  to  discuss  the  Riemann  Problem. 

The  formula  for  the  solution  to  the  Riemann  Problem  is  different 
for  different  orderings  of  the  initial  data  constants  u^,  uy  u^,  and  u^. 
This  indicates  that  twenty-four  cases  must  be  considered;  however  reflec¬ 
tions,  inversions,  and  reflected  inversions  of  cases  previously  discussed 
will  only  be  indicated  as  such,  and  will  not  be  discussed  in  detail.  This 
reduces  the  number  of  cases  to  be  discussed,  to  eight.  The  formula  to  be 
given  consists  of  a  picture  for  each  case,  with  labels  such  as  RXCu^.u^J, 
SYCu^.u^l,  and  rCu^u^l  given  to  parts  of  the  picture.  The  interested 
reader  may  then  refer  to  the  formula  given  earlier  in  this  chapter  for 
each  of  these  phenomena. 

Case  1.  Uj  <  u^  <  u^  <  u^.  At  t  =  1  the  solution  looks  like  Figure  4. 
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Figure  4.  The  Case  u,  <  u~  <  u,  <  u, 

i  2  3  A 

In  this  case  we  have  two  shock  waves  extending  from  the  intersection  of 
SY[u^,u^J  and  SXfu^jU^J  to  the  point  P,  and  from  the  line  x  =  f'(ii^)t, 
y  =  y(u^)t  to  P.  The  following  theorem  concerning  these  shock  waves  is 
due  to  Guckenheimer  [5]. 


Theorem.  A  shock  surface  between  two  constant  states  a,  b,  in  a  solution 
to  the  Riemann  Problem  lies  in  a  plane  which  passes  through  the  line 


Thus  the  two  shock  waves  extending  to  P  are  planar.  We  shall  call  these 
shock  waves  "4,-shocks." 

Now  from  (2.15)  we  deduce  that  P  has  coordinates 


Note  that  since 


u^,  we  have 


f(u2)  -  f(U[) 
ll2  "  U1 


f (u^)  -  f(Uj) 


<  f(0; 


(2.16) 


these  are  the  x-coord inates  of  SXlu  ,u,  ),  P,  and  the  left  edge  of  RX[u  ,u,  i, 

2  1  3  4 

respectively.  Thus  both  f-shocks  may  be  parameterized  by  s  and  t  with  the 
equations  x  =  f'(s)t  and  y  =  iKs)t.  In  case  u^  =  or  u2  =  then  either 
the  left  or  right  H'-shock,  respectively,  is  not  present,  because  both  end 
points  of  that  particular  shock  are  identical. 

Using  the  method  by  which  (2.12)  was  derived,  one  may  derive  the 
following  differential  equation  for 


<l>’(s)  =  f"(s) 


g(Uj)  -  g(u3)  -  lKs)(Uj  -  «3) 
■f(Ul)  -  f(u3)  -  f,(s)(ul  -  u3) 


(2.17) 


Note,  from  Figure  4,  that  the  right  'f-shock  meets  the  left  endpoint  of  the 
ft u4 * uf :  u3^  Shock  wave  continuously,  so  that  i{>(u3)  *  Comparing 

(2.17)  with  (2.12),  we  see  that  ^’(u^)  =  y’ (u.j) .  Thus  the  T-shock  meets 
the  Y-shock  with  first-order  smoothness.  We  may  also  deduce  that  for  any 
point  (f ' (v),  y(v) )  on  the  T[u^,Uj  ]  shock  at  t  =  1,  the  tangent  line  to 
the  curve  x  =  f'(s),  y  =  y(s)  passes  through  the  point  P  given  by 


f(v)  -  f(uJ)  g(v)  -  g(Uj ) 


(2.18) 


One  may  also  deduce  this  by  rewriting  (2.12); 


g(uL)  -  g(s) 


=  y'(s)  _ 
dx  f"(8) 


U1  -  S 


f^)  -  f(s) 

_______ 


(2. ly) 


f'(s) 
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Case  2.  u.  <  u.  <  <  ti_.  This  is  the  reflection  across  the  line  y  =  x 

-  1  4  J  2 

of  Case  1. 

Case  3-  u„  >  u.  >  u,  >  u„.  This  is  the  inversion  of  Case  1. 

-  .3  4  1  2 

Case  4.  u.  >  u„  >  u,  >  u. .  This  is  the  reflected  inversion  of  Case  1. 
-  3  2  14 

Case  5.  <  u^  <  <  u^.  At  t  =  1  the  solution  looks  like  Figure  5. 

Note  that  the  right  edge  of  RX[u2»u^]  has  equation  x  =  f'(u^)t. 
Since  f"  >  0  and  f’(u^)  <  f'(u^);  furthermore,  x  =  f'(u^)t  is  the 

equation  for  the  left  edge  of  RX[u^,u^J.  Thus  RX[u2>u^j  lies  completely 
to  the  left  of  RX[u^,u^j  in  this  case. 

The  point  P  has  the  same  coordinates  as  in  Case  1. 


1  5 


Case  6.  u,  <  u,  <  u.  <  u„.  This  is  the  reflection  of  Cose  5.  Note  that 
-  4  13  2 

Cases  5  and  6  are  invariant  under  inversion. 

Case  7.  <  u^  <  <  u^.  At  t  =  1  the  solution  looks  like  Figure  6. 

In  this  case,  since  u.^  <  u^,  and  f"  >  0,  we  have  that  f'(u^)  •  f'(Uj) 

hence  RX[u„,u,  1  overlaps  RYlu„,u(  I.  Also  the  shock  waves  rTu,,u,J  and 
2  1  3  4  4  1 

FlCu^.u^J  are  not  truncated  in  this  case.  The  point  has  coordinates 
(£ ' (Uj) ,g' (Uj)) ,  and  ().,  =  (f’Cu^),  g'fu^)). 

Case  8.  u.  <  u„  <  u.  <  u„.  This  is  the  reflection  of  Case  7.  Note  that 
4  3  1  2 

Cases  7  and  8  are  invariant  under  inversion. 


Figure  b. 


The  case  u„  < 


<  Uj_  < 


“I  < 


«3  <  u2  <  u4. 


Case  9. 


At  t  =  1  the  solution  looks  like  Figure  7. 


Figure  7.  The  case  <  u^  < 

The  point  P  has  the  same  coordinates  as  before. 

Case  10.  u.  <  u„  <  u,  <  u„.  This  is  the  reflection  of  Case  9. 

-  13  4  2 

Case  11.  u„  >  u,  >  u,  >  u„.  This  is  the  inversion  of  Case  9. 

-  3  14  2 

Case  12.  u„  >  u.  >  u„  >  u. .  This  is  the  reflected  inversion  of  Case  9 
3  12  4 

Case  13.  <  u^.  At  t  =  1  the  solution  looks  like  Figure  8 

Case  14.  <  u4  <  U2  <  l,3'  This  is  the  reflection,  and  also  the 


inversion,  of  Case  13. 


SX[u2,u1] 


SY[u3,u2] 


SY[u4,Ui] 


SXfuj.u^] 


Figure  8.  The  case  u^  <  u2  <  <  u^ 


Case  15.  u_  <  u  <  u,  <  u, .  At  t  =  1  the  solution  looks  like  Figure  9. 
-  I  3  A  1 

Here  Q  has  coordinates  (f'(u3),  g'Cu^)).  Between  the  points  Cj(f'(u4>, 
g'(uA))  and  C2  =  (f'(Uj),  g'(u^)),  the  rarefaction  waves  RXru^.u^]  and 
RYtu^.u^  I  meet  along  the  surface  A,  which  may  be  described  by  the  equations 
x  =  f’(s)t  and  y  =  g'(s)t,  for  u^  <  s  <  u^  and  t  >  0.  The  plane  sections 
where  u  =  s  in  each  wave  meet  along  the  line  x  *  f'(s)t,  y  =  g'(s)t. 

Thus  the  solution  is  continuous,  though  not  differentiable,  along  this 
surface.  One  may  check  that  the  entropy  condition  is  satisfied  near  A. 


Case  16.  u,  <  u„  <  u„  <  u. .  This  is  the  reflection  of  Case  15. 
- 4  3  2  1 


Case  17.  u,  >  u^  >  u2  >  u^. 


This  is  the  inversion  of  Case  15. 


RX[u2,Ul] 


lc2  RY[u4,Ul] 


SY[u3,u2] 


ri[u2>u3] 


RX[u3\u4] 

Figure  9 .  The  case  u2  <  u3  < 


Case  u„  >  u  >  u.  >  u_.  This  the  reflected  inversion  of  Case  15. 

-  2  14  3 

Case  J  9.  u„  <  u,  <  u,  <  u..  At  t  =  1  the  solution  looks  like  Figure  10. 

- , —  Z  4  J  1 

Case  20.  u^  <  u2  <  <  u^.  This  is  the  reflection  of  Case  19. 

Case  21.  u.  >  u„  >  u,  >  u_.  This  is  the  inversion  of  Case  19. 

-  4  2  13 

Case  22.  u„  >  u.  >  u.  >  u„.  This  is  the  reflection  inversion  of  Cise  19 
-  z  4  1  3 

Case  23.  u3  <  u2  <  u^  <  u^.  At  t  =  1  the  solution  looks  like  Figure  11. 
Here  has  coordinates  (f'(u3),  g'(u3)).  Also  Q2  =  (f'(u2),  g'(u0)), 

Q3  m  g’(u4)),  q4  =  (f’Cu^,  g'Cuj)). 

Case  24.  u_  <  u,  <  u„  <  u..  This  is  the  reflection  of  Case  23.  Note 
-  3  4  2  1 

that  Case  23  and  24  are  invariant  under  inversion.  Also  note  that  Cases 


23  and  24  are  the  only  cases  with  continuous  solutions. 


3.  Verification  o f  the  e nt ropy  conditio n . 

We  will  now  prove  three  theorems  showing  that  under  certain  condi¬ 
tions,  our  solutions  are  single  valued  and  satisfy  the  entropy  condition. 
Theorem  1  treats  the  case  where  f  =  g  and  is  convex.  It  is  easily  seen 
that  this  also  includes  the  case  f"  =  eg".  Theorems  2  and  3  concern 
pertubations  away  from  this  case. 

Theorem  1 .  If  f  C-  g,  then  the  constructions  in  §2  for  Cases  1-24 
define  functions,  and  these  functions  satisfy  the  entropy  condition. 

Remark.  In  Cases  9-12  and  19-22  it  is  conceivable  that  our  solution 
may  fail  to  be  single  valued  due  to  overlapping  T  shocks,  as  illustrated  in 
Figure  12.  Therefore  it  is  neccessary  to  prove  that  this  does  not  occur. 


y 


l 
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Proof .  We  need  to  prove  chat  Che  entropy  condition  is  satisfied 
for  SX,  SY,  T,  and  'P.  However  SX  and  SY  are  one-dimensional  shocks 
and  the  entropy  condition  is  known  to  hold  for  them.  The  proof  for 
r[v,  w]  is  as  follows:  To  verify  the  entropy  condition  we  require  n, 
the  normal  vector  to  r[v,  w],  to  be  oriented  towards  the  side  of  the 
shock  surface  where  u  is  larger.  In  this  situation  this  means  that 
n  must  "point  in"  to  the  rarefaction  wave.  Since  the  rarefaction  wave 
lies  in  the  region  y  <  y(s)t,  we  must  have  n^  <_  0.  Since  f"  >  0,  this 
means  we  should  choose: 

n  -  (f"(s) y(s)  -  f ' (s)y' (s) ,  y’(s),  -  f"(s)).  (3.1) 

Thus,  we  must  verify 

(k-w)(f"(s)  Y(s)  -f(s)  Y'(s))  +  (f(k)-f(w))(Y'(s)-f"(s))  ?0  (3.2) 

for  all  w  <  k  <  s  <  v 

Lemma  1 . 1  When  f  =  g,  y'(s)  <  f"(s)  for  s  >  w. 

Proof.  Using  (2.13),  and  using  the  fact  that  f"  >0,  we  have 
that  y(v)  <  f'(v).  Furthermore  y  =  f '  is  a  solution  of  (2.12).  The 
uniqueness  of  solutions  for  s  >  w  implies  that  y(s)  <  f'(s)  for  all 
s  >  w.  Using  (2.12)  and  (2.13),  one  notes  that  y'(v)  -  0  <  f"(v) ,  and 
also  that  y'(s)  =  f"(s)  implies  f(w)  -  f(s)  -  y(s)(w-s)  *  f (w)  -  f(s)  - 
f'(s)(w-s);  this  in  turn  implies  that  f'(s)  ■  y(s)  when  s  f  w.  Since 
Y(s)  <  f'(s),  we  must  have  that  y'(s)  <  f"(s)  for  s  >  w. 

We  now  verify  (3.2).  If  the  left  hand  side  of  (3.2)  is  consid¬ 
ered  as  a  function  F(k),  then  one  checks  using  Lemma  1.1  that  F"(k)  >  0; 


22 


also  F(w)  »  0.  Furthermore,  F(s)  *  0  since  in  (3.2)  one  may  replace 
all  w's  by  s's,  using  (2.1).  Therefore  F(k)  >  0  for  w  <  k  <  s  <  v; 
this  is  (3.2). 

The  entropy  condition  for  reflected  and  inverted  T[v,  w]  shock 
waves  may  similarly  be  verified. 

It  remains  to  verify  the  entropy  condition  for  Y-shocks.  These 
appear  in  two  different  contexts: 

(1)  Tangential  shocks,  that  is,  ¥  shocks  which  are  tangential 
to  T  shocks.  Such  a  V  shock  has  the  same  normal  vector  n  as  does  the 
T  shock  at  the  point  of  tangency,  and  also  the  same  values  of  u  on 
either  side  of  the  shock.  Therefore  the  entropy  condition  for  a  tang¬ 
ential  ¥-shock  is  equivalent  to  the  entropy  condition  for  the  corres¬ 
ponding  T  shock  at  the  point  of  tangency. 

(2)  Non-tangential  shocks.  One  may  check,  using  the  ordering 
of  u^,  u2>  u^»  and  u^ ,  and  the  convexity  of  f  and  g,  that  the  inter¬ 
section  of  any  non-tangential  ¥-shock  surface  with  the  plane  t  =  1 

is  a  line  segment  with  negative  slope,  as  depicted  in  Figures  4  and 

8  •  Furthermore  any  non-tangential  ¥-shock  is  a  shock  between  u^  and 

u^,  with  u^  >  u^,  with  the  region  u  =  u^  above,  and  to  the  right  of 

the  shock.  Thus  n  ■  (n  ,  n  ,  n  )  with  n  and  n  negative.  The 

t  x  y  x  y 

entropy  condition  becomes 

(k  -  u3)nt  +  (fOO  -  f(u3))(nx  +  ny)  >_  0,  (3.3) 

for  u,  >  k  >  u  .  Since  n  +  n  <0,  and  f  is  convex,  the  left  side  of 
3  1  x  y 

this  inequality  has  negative  second  derivative  with  respect  to  k,  and 
by  (2.1),  is  zero  when  k  ■  u^  or  k  «  u  . 


Hence  it  is  positive  for 


r 
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values  of  k  between  and  u.^. 

To  verify  that  our  solutions  are  single  valued  it  suffices  to  show 
that  overlap  of  F-shocks  does  not  occur.  In  the  proof  of  Lemma  1 . 1  we 
saw  that  y(s)  <  f'(s);  thus  both  T  shocks  lie  on  opposite  sides  of  the 
curve  x  =  "’(s),  y  =  f'(s). 

Theorem  2.  Let  u^,  ,  u^,  u^  he  such  that  our  proposed  solution 

to  (1.1),  (1.2)  contains  no  complete  F,  FR,  FI,  or  TIR  shock,  that  is, 

let  us  consider  only  Cases  1-6,  9-14,  23,  and  24  of  §2  .  Let 

M  =  max  u.,  and  m  =  min  u..  Then  for  any  given  function  h  such  that 
l<i<4  1  l<i<4  1 

h"  >  0,  there  exists  c  >  0  such  that  whenever  [j  f  -h||  „  <  e,  and 

C  [m,M] 

| [ g  —  h | [  <2,  our  construction  for  the  solution  to  (1.1),  (1.2)  in 

C  [m,M] 

these  cases  defines  a  function,  and  this  function  satisfies  the  entropy 
condition. 

Proof.  As  noted  before,  it  suffices  to  prove  that  in  our  solution 
to  the  perturbed  equation,  the  entropy  condition  holds  for  truncated 
T-shocks,  and  for  non-tangential  ^-shocks. 

Recall  that  the  entropy  condition  states  that  for  T[v,w:  p]  we 
must  have: 

(k  -  w)n  +  (f (k)  -  f (w) ) n  +  (g(k)  -  g(w))n  >  0,  (3.4) 

t  x  y 

or 

W(k, s)  =  (f'(s)y(s)  -  f'(s)y'(s))  (k  -  w)  (3.5) 


+  y'(s)(f(k)  -  f(w))  -  f"(s) (g(k)  -  g(w))  >  0, 
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for  all  w  <  k  <  s  S  v,  and  w  •-  p  ■  s.  The  function  y  satisfie  the  initial 

value  problem  (2.12),  (2.13).  Since  f  is  convex  (for  £  sufficiently  small) 

and  s  ^  p  >  w,  the  denominators  in  (2.14),  the  formula  for  y,  are  bounded 

away  from  zero,  and  the  bound  depends  only  on  h,  c,  v,  w,  and  p.  Thus  we 

2 

may  conclude  that  the  map  (f,g)  -*■  /  is  continuous  from  {fjf  £  C  [p,vl, 
f"  >  0}  x  fg|g  £  C2[p,vJ,  g"  >  0}  to  C*[p,v.1. 

Next  note  that 

=  y'(s)f"(k)  -  f"(s)g"(k) ,  (3.6) 

3k2 


and  note  that  W(s,s)  =  W(w,s)  =  0,  by  (2.1).  When  f  3  g  =  h  we  know  by 

32W 

Lemma  1.1  that  y'(s)  <  f"(s)  for  p  <  s  <  v.  Thus  — -  <  0  for  w  <  k  £  s  <  v 

w  <  p  <  s,  when  f  =  g  =  h.  Moreover  -2  depends  continuously  on  f  and  g 
2 

in  the  C  topology  on  f  and  g.  Therefore  there  exists  e  >  0  such  that 

32W 

fff-h||  ?  <  e  and  ||g-h||  9  <  £  imply  — ~  >  0,  and  it  follows  that 
C  C  3k 


r[v,w:  pj  satisfies  the  entropy  condition. 

Since  y(s)  <  f'(s)  when  f  =  g  =  h,  for  £  sufficiently  small  we  must 
have  y(s)  <  g'(s),  thus  in  Case  9-12  the  two  i'- shocks  x*'  on  opposite 
sides  of  the  surface  x  =  f'(s)t,  y  =  g’(s)t.  Thus  overlap  does  not  occur 
for  e  sufficiently  small. 

Finally,  for  non-tangential  'f-shocks  we  must  show: 


(k  -  u3)nt  +  (f(k)  -  f(u3))n^  +  (g(k)  -  g(u3))ny  5  0  (3.7) 

for  u.  5  k  <,  u_.  Since  n  <0  and  n  <  0  as  we  saw  in  the  previous  sec- 
13  x  y 

tion,  and  since  f"  >  0  and  g"  >  0  for  £  sufficiently  small,  the  left  hand 
side  of  (3.7)  has  negative  second  derivative  with  respect  to  k,  and  equals 
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zero  when  k  =  Uj  or  k  =  u^.  Thus  (3.7)  holds  for  values  of  k  between  u^ 
and  u^. 

Theorem  3.  Suppose  u^,  u^,  u^,  and  are  such  that  our  proposed 

formula  for  the  solution  to  (1.1),  (1.2),  contains  some  complete  r[v,w], 

raCv.w),  ri[v,w],  or  riRfv.w]  shock,  that  is,  let  us  consider  Cases  7,  8 

and  15-22;  suppose  also  that  f"(w)  =  g"(w)  and  f'"(w)  =  g" ' (w)  at  all  w 

which  occur  as  above.  (Note:  w  is  always  either  u^  or  u^).  Let  m  and 

M  be  as  defined  in  Theorem  2.  Then  for  any  given  function  h  such  that 

h"  >  0  and  h"(w)  =  f"(w),  h'"  (w)  =  f"'(w),  there  exists  e  >  0  such  that 

if  ||f-  hjj  .  <  E,  ! j  s  —  h j  1  ,  <  e,  then  our  proposed  solution  is 

C  fm,Mj  C  [m,M] 

single  valued,  and  it  satisfies  the  entropy  condition. 


Proof .  The  proof  consists  of  several  lemmas. 

Lemma  3.1.  The  map  T:  C^[v,vJ  *  C4[v,w]  ■+  Cl[v,v],  T(f,g)  =  y, 
where  y  satisfies  (2.12),  (2.13),  is  continuous  for  f"  >  0,  g"  >  0.  The 

J 

map  S  (f,g)  =  -  T(f,g)(s)  =  yf(s)  satisfies 

&  Lie* 


UDSs)(f  <-  c\s  -  u\  li  (p,  q)  II  , 

1  C  [v,w] 

for  all  p,  q  in  C  *[v,w.J  .such  that  p"(w)  =  q"(w)  *  p"'(w) 
where  C  >  0  does  not  depend  on  p,  q,  or  s. 


(3.8) 

q" '  (w)  =  o, 


Proof .  From  (2. 1  A)  we  have 
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=  (def.)  T  (f,g).  Clearly  Is  a  continuous  and  differentiable  mapping 
from  {f  j  f  £  C4rv,wj,  f"  >  0)  x  {gjg  £  CA[v,w],  g"  >  0}  to  JR,  for  s  £  (v,v 
However  it  is  necessary  to  show  that  ^  is  bounded  independent  of  s, 

and  independent  of  the  choice  of  (f,g)  from  a  neighborhood  of  (h,h). 

Define  Tw(f,g)  to  equal  g'(w).  We  prove  that  s  Tg(f,g)  is  con¬ 
tinuous  at  w,  as  follows.  We  have: 


Y(s)  -  T  (f.g)  -  ~  6(W)- 

C  r-  r  t 


(3.10) 


(r  x-umw-o  y dr 


S  /*  f"(0)(0-w)  d0. 


(w-r)  ‘ 


Note  that 


r  fH(z)(w-z)  dz 
a  f*  f  "(9)  (0-w)  d0< 


(3.  11) 


At  2  f*  £" '  (  0)  ( 9-w) 1  d0 

exp(  7^1  +  - 2 - 7z - 2 - 

Vs  (w-z)(f"(z)(z-wr  -  l  f'"(0)(0-wr  d 0) 

w 


<s-w>2  a  /;  A 

- -  exp  r  -  Hz  . 

(r-w)  \s  (w-z)  /z  f"(0)(9-w)  d0  / 

w 


|y(s) 


;(s)  -  g(w) 


(3.12) 


<  (s-w)' 


. » exp  0 


2Hf,"L  V  II*" II . 
a  3(inf(f,,))  7  (-w)2 


(s-w)z  ( 


2  r  l 


s-w  V 


!l8"L  ^llf-IL(v-w)’ 

J  2  exp  \  3  ( inf (f ") ) 


r 
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and  so  we  see  that  lim  Y(s>  -  g'(w);  thus  s  -  T  (f,g)  is  continuous  on 


S~*W 


l"v,wj.  The  sup  norms 
over  I v, w  I . 

Next , 


5  and  infimums  used  above  and  henceforth  are  all  taken 


,  ,  s  q(s)  -  q(w) 

(DTs)  (F,g)(p’ll>  "  8  -  » 


(3.13) 


exp 


f"(2)(w-2) 

s  !Z  f"(0) (0-w)  d0 
*  w 


dz 


r 


p"(z)(w-z)  JZ  f"(0)(0-w)  d0  -  f"(z)(z-w)  JZ  p"(6) (0-w)  d0 


dz 


/Z  f"(d)(Q~w)  dd 


r  de  +  fr 


w  (r-w)' 


d9 


w  (r-w)‘ 


dr 


Thus 


I  (DTs)(f>g)(p,q)  I  <  llq'l 


(3.14) 


/•V  ,  s2 

(s-w) 


s  (r-w) 


/2||£m,il  (v-w>n 

2  exp  \  3( inf ( f ")) 


(Jr  tj2  f " (6) ( 0-w)  d9 ] ~L  [(p"(w)  +  JZ  p"'(0)d0)(w-z)  JZ  f"(6) (O-w)dO 
s  w  w  w 

VL 


-  ~  f"(z)(z-w)(p"(w)(z-w)“  -  /Z  p"’(0)((0-w)z  -  (z-w)Z)dO)]  dz)  - 

L  4  W  ? 

l  q"(w) (r-w) 2  -  II  q" ' (9) [ (0-w) 2  -  (r-w)2]  d0 


(r-w) ‘ 


dr. 
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Since  only  pc  r  t  n  i'l>.  i  L  i  ons  fixing  f"(v)  ,  g"(v), 
are  considered,  we  have  p"(w)  ~  i|"(w)  ;  p"'(w)  -  n"' (v)  =  . 


(3.17)  !  (DTs)  (f  }  (p.<|)  I  ■  !|q'!l 


rj  ,  .2 

,  (s-w) 

+  - exP 

Js  (r-wV 


'2;|r,’||„,  (v-w)' 

"ITinYfr-T)" 


('i ♦  ■!) 

,'!c.  !”■+  i!q»'  ii  (r_w) 

( inf (f ") ) 2 

v  i  ) 

2  +  '■  *  3 

dr 


.1  111'  II  +  exp  (- 


(v-w)' 
3(  iiil'Tf  ")  j 


2  ^3'  '-s-w  v-w' 


( Inf ( r") ) 1 


+  ^  th"  ’  II  On  (v-w)  -  1  n ( s-w) ) 


J 


(s-w) 


nnd  so  wo  see  that  ||(l)T  ).  .  |(  is  independent  of  s,  w  <  s  <  v,  and 

s  1 1 .  g)  _  _ 

also  Independent  of  the  choice  of  (f,g)  from  a  sufficiently  small 

neighborhood  of  (h,h) .  Therefore  by  the  Mean  Value  Theorem  the  map 

(f.g.s)  >  y(k)  is  locally  l.ipschit/.  for  each  S,  with  a  l.ipschitz  con 

slant  independent  of  s.  Thus  (f>g)  ►  v  is  a  continuous  mapping  from 


3  3 

C  [  w,v]  x  C  lw,v]  C[w,v|. 


Next , 


(3.18)  s  ( f , g)  —  y * ( ; . )  p’ r. 

(s-w)2 


g(y) .at?)  .“  c’ (s).(w-s)_  _  fs  j  rr _ f "(?.) (w-z) 


(s-w) ' 


Jv  Vs  [*  ("(C)  (H-w)  do. 


-£"(,)(»-»)  \  II  <» 

I  >  dl  . 

f"O)(0-w)  do/  (r-w)' 
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Differentiating  with  respect  to  f  and  g,  we  have 


<DSs,(f,E)(p'<l) 


(3.17) 


(w-s) 


[p"(s)  Is  f"(0)(b-w)  d9  -  f"(s)  /*  p"(6) (0-w)  de ] 
w  4  w 


I/®  f "(e) (6-w)  der 


exp 


-  dz)  <1  *"(9>ilirl ”  Jt 


s  Jz  f"(0) (e-w)  d 0 


( r-w)  ‘ 


f(s)(w-3) 

Is  f "(0) (0-w)  d0 


cr 


exp 


f**(z)  (w-z) 

[Z 

1  w 


s  Jz  f"(0) (0-w)  d0 


dz 


fz 


f "(0) (0-w)  d0  "  (w-z) 


V^w 


9  ~2  (w-z)  |V(z)  f" 


1  ( 0 ) ( 0— w)  d0 


-  f"(z) 


p"(0)(0-w)  del  d^  fr  da  +  f 


dr 


;w  (r-w)“  ;w  (r-w) 

(w-s)  [f3  p,"(0)dfl  /s  f"(0)  (0-w)d0  +  /*  p"  ’  (0)((0-w) 2  -  (s-w)2)d9] 

W  W  w _ _ _ _ 


f/®  f "(0) (e-w)  d0] 


/  s-w-i  2 

1 - J  exp 

'-r-w-' 


fz  f " '  ( 0)  ( 9— w)  2  d0 
w _ _ 

s  (w-z)  /z  f"(0) (0-w)  d0 


dz) 


r  g"(0) (0-w)  d9 
w  (r-w) 


dr 


+  f"(s) (w-s)  [S 

f3  f"(0)(O-w)  d0  Jv 


exp 


r  /Z  f " ' (9) ( 0-w) 2  d0  \  , 

- -  dz  (*=«)  2 

s  (w-z)  JZ  f"(0) (0-w)  d9  y  r  W 


r  r  (z 

(w-z)  f"(9)(0-w)  d0 

S  l  'W 


•]- 


Jd 
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p"'(9)  dd 


w 


f "(9) (6-w)  d6 


w 


+  f"(z)  P  P'^—  ( (0-w) 2  -  (z-v)2)  d9  j  dz 
■*w 


r  g"(9)  (6-w)  de  _  [r  q '  "(9X(9-w) 2  -  (r-w)2) 


w  (r-w)‘ 


w 


2 (r-w)' 


'J' 


d9  dr. 


Note  that,  since  p"'(w)  =  0, 


p"'(0)  d0  =  - 


pM" (0) (( 0-w)  -  (s-w))  dO, 


(3.18) 


and 


f 

Jw 


pM 1 (0) (8-w) L  d0 


■D 


p”"(e)  ((z-w)3  -  (9-w)3)  de. 


(3.19) 


Thus 


<DS3)(£,g)(p.q) 


[<-,[ 


rs 


f "(9) (9-w)  d0 


w 


0 


-2 


(3.20) 


p ""(0)  ((s-w)  -  (9-w))  d0 


f"(9) (0-w)  d0 


w 


f  P,,M(0) 

r  3 

(s-w) J 

6 

<e-„)3  _  (s-v)2  ((s.v)  .  (e.u))l 

de 

'W 

•J 

— - 

(S)2  «p 


r  j;  r(W-.)!  de  ^  r  8..(e)(e-M) 

s  (w-z)  /*  f"(0) (9-w)  d0 


dz)  •  ft-vr/Ar  d0  dr 
w  (r-w)Z 


V 
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Thus  j|  (DS  )  (f,g)j|  <  C  is  -  w|  for  some  C  >  0,  where  C  depends  .l'y  on 
i|fM,IL  •  inm"),  i|g"!L  •  ;ind  (v,w).// 

Lemma  3.2.  In  a  I’i  v,w j  shock,  in  the  case  where  f  =  g,  Y"(w)  <  f"'  (  ■ 


Proof.  Recall  that  when  f  =  g,  f'  Is  a  solution  to  (2.12), 


Thus 


y(s)  =  f'(s)  +  a(s),  where  a  satisfies 

f”(s)(s-w) 


a' (s)  -  a(s) 


(3.22) 


fS  f"(0)(e-w)  d9 

*  W 


Thus 


(s-w)  ~r~  In  (a(s))  = 
as 


f"(s)(s-w)^ 

JJ  fM(0)<e-w)  de 


(3.23) 


Note  that 


lim 

s-*-w  fS  f"(e)(0-w)  d0 

*  u 


(3.24) 


Thus  one  may  write 


o(s) 


C ( s-w) 2 


exp 


f"'(0)(e-w)2  d0  \ 

- £ -  dz 

(z-w)  fl  f"(9)(0-w)  d0  / 
w 


(3.23) 


One  may  observe  that  o  has  2  continuous  derivatives  at  w  if  f  has  3  con¬ 
tinuous  derivatives,  and,  0"(w)  <  0  if  and  only  if  C  <  0. 


Thus  Y"(w)  <  f"  (w)  if  and  only  if  "V  <  f  ’ .  However 


n 


y/v\  -  ISyJi  ~  <  f'tv)  since  f"  <  0. 

We  can  now  prove  that  for  sufficiently  small  £,  .If-*',,  ,  <  t  and 

H  g  —  h 1 1  <  e,  together  with  the  hypotiiesos  of  Theorem  3,  implv  that  over- 

C4 

Lap  does  not  occur.  It  suffices  to  show  that  for  t  sufficiently  small, 

Y(s)  <  f ' (s)  ,  w  <  s  <  v.  From  this  it  follows  that  the  line  y  =  x  separates 

the  two  T-shocks  in  Cases  19  through  22. 

Since  when  f  =  g  =  h,  y(s)  <  f'(s)  for  w  <  s  <  ,  and  since 

y  depends  continuously  on  f  and  g,  we  have  that  for  any  6  y  0  there 

exists  e  >  0  such  that  ||f — h ||  ,  <  e  and  ||g-h||  ,  <  e  Imply  y(s)  <  f  ts) 

C  C 

for  w  +  6  <_  s  <_  v.  For  s  near  w,  and  f  =  g  5  h, 

f"(s)  -  y ' (s)  -  (fM,(w)  -  y"{w) ) ( s-w)  +  o(s-w) .  (3.2b) 

Since  by  Lemma  3.2  f'"(w)  -  y"(w)  >  0,  for  6  sufficiently  small  we  have 
f"(s)  -  y ' (s)  >  CQ(s-w)  (3.27) 

for  some  Cq  >  0  and  for  w  <_  s  6 .  By  Lemma  3.1  and  the  Mean  Value 

theorem,  for  ||f-h||  ,  <  e  and  ||g-h||  ,  <  e 
C  C 

f"(s)  -  y*  (s)  >•  (CQ  -  eC)  (s-w)  (3- 23) 

for  w  <_  s  <  6.  For  e  sufficiently  small  Cq  -  eC  >  0,  and  hence  f"(s) 

>  y'(s)  for  w  <  s  £  6.  Since  y(w)  =  f’(w),  the  result  follows. 

To  finish  the  proof  of  Theorem  3,  it  must  be  shown  that 

W(k,s)  =  (f"(s)y (s)  -  f ' (s)y ' (s) )  (k-w)  (3-29) 

+  y ' (s) ( f Ck)  -  f(w))  -  f"(s)(g(k)  -  g(w) )  >0, 


Recall  that  W(s,s)  =  W(w,s)  =  0,  and. 


for  all  k,  s,  w  <.  k  <.  s  <  v. 

32W 

when  f  =  g,  — j  *  f"(k) (y ' (s)  -  f"(s))  which  is  less  than  zero  for 

3k 

w  <  s  <  v.  In  fact,  for  s  close  to  w, 


f "(w) (y"(w) 


-  f'"(w))(s-w)  +  o(s-w). 


O.30) 


On  the  other  hand. 


°(f,g) 

'  °(f,g) 


(Y’)(p,q)f”(k)  +  y' (s)p"(k)  -  p"(s)g"(k)  -  f"(s)q"(k) 


Thus 


(3.32) 


<  O(s-w)  ||(p,q)  l|  ,  f”(k)  +  |  Y ' (s)  P"(k)  -  p"(s)g"(k)  -  g"(s)q"(k)| 

C 

Since  we  are  considering  only  those  tangent  vectors  p  and  q  such  that 
p"(w)  ■  q" (w)  -  0,  we  have  that  |p"(k)|,  jp"(s)|,  and  |q"(k)|  are 
less  than  ^||p||  ^  +  ||qj|  |  s-w|  .  Thus  for  some  >  0, 


|D(f,g)  7T)(p,q)i  -ci  *  |s"w!  *  II 


Ok 


(3.33) 


2 

aw 


Thus  for  f  =  g  5  h,  — j  s  Cq  •  |  s-w|  for  some  >  0  ;  and  for 

3  le 
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f  =  h  +  ^p,  g  =  h  +  £q,  where  p  nnd  q  satisfy  p"(w)  =  q"(w)  =  p'"  (w)  = 
q"' (w)  =  0, 

-h  *  (E(C  +  C ,)  ||  <P,q)  |t  4  -  C  )|s  -  w|  <0  (3.34) 

1  C  u 

for  sufficiently  small  £.  Since  under  these  conditions  W  is  convex  down, 
and  W  =  0  at  k  =  s  and  k  =  w,  one  concludes  that  W  ^  0  for  w  <  k  <  s  <  v  ./// 

Remark.  In  Theorems  2  and  3  one  may,  of  course,  fix  f  and  perturb 
only  g,  or  the  other  way  around.  In  this  case  one  chooses  h  =  f,  or  h  =  g, 
respectively. 

4 .  A  Counterexample. 

OO 

The  following  is  an  example  of  a  C  one  parameter  family  (f,ft)  of 

OO 

pairs  of  C  functions  such  that  f  =  and  such  that  for  certain  initial 

data,  the  solution  given  in  §2  does  not  satisfy  the  entropy  condition  for 

l-  >  0.  However,  we  also  give  the  correct  entropy  solution  for  this  example. 
2  2  3 

Let  f(s)  =  s  ,  and  f£(s)  =  s  +  es  .  Consider  the  initial  value 

problem 

3  3  3 

u(t,x,y)  +  —  f (u(t ,x,y) )  +  fe(u(t,x,y))  =  0  ,  (4.1) 

u(0,x,y)  =  6  for  x  >  0  and  y  <  0  , 

=  0  otherwise, 

where  6  >  0.  Note  that  f"  >  0  on  TO, 61  for  e  >  -1/(36).  The  solution 
given  in  f)2  is  described  in  I’igure  13. 

2  2  3 

Proposition.  If  f(s)  =  s  ,  g(s)  =  f£(s)  =  s  +  ts  ,  then  for  no 
choice  of  t  >  0,  and  6  such  that  0  <  6  does  the  shock  f[6,0]  satisfy  the 


entropy  condition,  near  the  line  x  =  0,  y  =  0. 


Also 


Y ' (s)  =  2  -  4es  ln(s)-2es  +  2s  (-  -J  +  eC21n(6)  +  1) ) .  (4.4) 

The  shock  wave  T[6,0]  satisfies  the  entropy  condition  if 

W(k,s)  =  (k-0) (2y(s)  -  2sy’(s))  +  (k2-0)y'(s)  (4.5) 

-  (k2  +  ek2)  •  2  >  0  for  0  <  k  <  s  <  6. 


,)W 

Note  that  —  (0,s)  =  2y(s)  -  2sy'(s).  Solving  y(s)  -  sy'(s)  =  0, 
we  have 


2s  -  2es“ln(s)  +  s~  (-  4  +  e(21n(a)  +  1)) 

a 


-■ 


2s  -  4es2ln(s)  -  2es2  +  2s2(-  4  +  s(21n(5) 


+  1))1-  0; 
J 


2es2ln(s)  +  2es2  -  s2  (  -  ^  +  e(21n(5)  +  1))  =  0; 


s2  (2eln(s)  +  2e  -  (-  4  +  e(21n(<5)  +  1))  =  0. 

c* 


(4.6) 


Thus  y(s)  =  sy'(s)  if  s  =0  or 


s  =  exp  (-  2^  (2e  -  (-  J  +  e(21n(5)  +  1))))  (4.7) 

-  (clef.)  Uq • 

Note  that  u  is  positive,  and  that  lim  u  =  0  for  0  <  6  .  Furthermore, 

U  C-+0  0 

sw 

for  0  <  s  <  Uq,  Y(s)  -  sY'  (s)  <  0,  hence,  —  (0,s)  «-  0  for  these  values 
of  s,  and  thus  W  <  0  for  these  values  of  s,  and  k  close  to  zero.  Thus  for 
any  choice  of  6  greater  than  0  our  solution  does  not  satisfy  the  entropy 
condition  for  small  £. 

To  keep  the  computations  relatively  simple,  we  give  the  correct 


entropy  solution  for  this  example  only  in  the  case  0=1.  In  this  case 
we  will  see  that  r[l,0:  uQ  I  satisfies  the  entropy  condition.  To  the  left 
of  x  =  f'(uQ)t,  a  new  rarefaction  wave  appears,  and  the  shock  wave  passes 
between  this  new  rarefaction  wave  and  RXfO.l].  At  t  =  1  the  solution  looks 
like  Figure  14. 


Let  us  call  the  new  rarefaction  wave  £,  and  the  continuation  of 
r[l,0:  Uq],  n.  Let  n  have  equations  x  =  f’(s)t,  y  =  w(s)t.  Below  fl,  the 
solution  u  equals  s  on  plane  sections  x  =  f'(s)t.  Above  fi,  u  =  v  on 
plane  sections 


r 


t<> 


y  -  f^(v)t 


w'(s) 

f"(s) 


(X 


t ' (v) t)  . 


(4.8) 


These  plane  sections  meet  the  shock  curve  at  x  *  f’(s)t,  y  *  w(s)t 
tangentially,  as  suggested  by  (4.8). 

We  now  prove  that  this  description  is  correct,  and  give  an  ex¬ 
plicit  expression  for  w  and  the  relationship  between  s  and  v. 

Since  SI  is  to  satisfy  Condition  2.1,  we  have 


w’(s)  -  f"(s) 


f  (v)  -  f.(s)  -  w( s) (v-s) 


f (v)  -  f(s)  -  f'(s)(v-s) 

Furthermore,  by  hypothesis,  u  =  v  along  a  plane  tangent  to  x  * 
y  =  w(s)t.  Thus  the  following  result,  due  to  Guckenheimer  [5] 
allow  us  to  get  a  different  expression  for  w'(s). 


(4.9) 

f ' (s) t, 
will 


Theorem.  If  u  is  a  solution  to  the  Riemann  Problem,  then  in¬ 
side  each  region  of  rarefaction,  the  surfaces  u  =  v  are  sections 
of  planes  passing  through  the  line  x  =  f'(v)t,  y  =  f^(v)t. 

From  this  result  we  get 


/ 


dx 


w' (s)  _ 
f"(s)  =  dx 


w(s)  -  f’(v) 
f'(s)  -  f ' (v) 


(4.10) 


Thus 


2  3  2  3 

V  +  -  s  -  £S  -  w(s)(v-s)  w(s)  -  2v  -  3e  v" 

2  2  .  .  .  2s  -  2v 

v-s  -  2s(v-s) 


w(s)  -  2v  -  3ev 


2  2  3  2  3 

(v  +  ev  -  s  -  es  -  w(s)(v-s)) 


(s-v) 


2(w(s)  -v-s  -  e(v^  +  vs  +  s^)). 


(4.11) 


AO 


And  thus 


2  2 
0  =  fcv  -  (2t-s)v  +  w(s)  -  2s  -  2ts  . 


(A. 12) 


So 


v  =  s  +  (3s2  -  (w(s)/c)  +  (2s/t_))1//2. 


(A. 13) 


Since  v  <  s  is  desired  we  choose  the  sign.  Now  by  (A.'O) 


.  ,  ,  w(s)  -  (2v  +  3tv  ) 

w’(s)  =  — - 

s  -  v 


(A.  14) 


„  ,  ,  ,  „  2  ,  2s  -  w(s)(l/2 

=  2  +  fits  -  4 e  (3s  +  - - - ) 


Now,  differentiating  (4.13)  with  respect  to  s,  and  substituting  (4.14), 
we  compute 


dv 

ds 


6s  + 


6s  + 


2  -w’(s) 


1  - 


2  (3s2  + 


(4.13) 


1  - 


1  .l( 


2  +  6es  -  4e 


/3s2  +  2s^_w(iT) 


-1 


2  /3s2  +  2s  -yisL 


Thus  v  =  -S  +  c.  Since,  as  depicted  in  Figure  14,  v  =  0  when  s  =  uQ, 
we  have  that  c  =  u^.  Thus  v  =  Uq  -  s,  and  we  expect  that  the  surfaces 
x  =  f’(s)t,  y  =  w(s) t ,  and  x  =  f'(v(s))t,  y  =  f^.(v(s))t  meet  when  s  =  v(s)  - 
Uq  -  s;  that  is,  at  s  =  Uq/2.. 

We  may  now  rewrite  (4.10): 


w'(s) 


w(s)  -  2v  -  3ev 


s  -  v 


(A.  In) 


w(s)  +  2s  -  2uq  -  3e(s  -  Uq)‘ 
2s  - 


Thus  we  have  a  linear  first  order  differential  equation  for  v.  The 


42 


we  have  that 


w(s) 


-t:(s 


uo/2)‘ 


+  (2  + 


3iuq)(s  - 


V2> 


uoa 


+  3iuq/4) .  (4 . 


19) 


Thus 


w'(s)  =  -2e(s - j)  +  2  +  3eUq,  (^ -  20) 

w"(s)  *  -2e  <  0,  for  e  >  0, 

and  we  see  that  the  shock  curve  y  *  w(s),  x  m  f'(s)  *  2s  is  concave 
down.  Furthermore  w'(ug/2)  =■  2  +  3eUg  **  2  +  6c(Ug/2)  =  f^(Ug/2). 

Thus  the  shock  curve  meets  the  curve  y  *  f^(v),  x  “  f'(v)  tangentially. 
To  verify  the  entropy  condition  for  ft,  it  is  necessary  to  show 

that 


W(k,s) *  (2w(s)  -  2sw' (s) )  (k-(ug-s) )  (4.21) 

+  (k2  -  (uQ-s)2)  w’(s)  -  (k2  +  Ek3  -  (uQ-s)2  -  e(:!q-s)3)2  ^  0, 

for  Uq-s  <  k  <  s,  (Uq/2)  £  s  Uq.  Substituting  (4.19),  we  have 

W(k,s)  =  (k  +  s  -  Ug) (2«s  -  2eUq)  +  (k2  -  Ug  +  2ugS  -  s2)[4eug  -  2tsJ 

r  3  3  2  27  (4.22) 

-  2[fk  -  eug  +  3tUgS  -  3£UgS^  +  esJ]  . 

It  is  now  easily  checked  that  W(s,s)  *  W(Ug  -  s,s)  ■  0;  this  verifies 
Condition  2.1,  the  Rankine-Hugoniot  condition,  for  ft.  Next 

(k,s)  ■  (2&s2  -  2euq)  +  2k  (4eug  -  2es)  -  6ck2 


(4.23) 


0  for  k=*--|s+-jUQ±-j(2s-  uQ) 


h  s  +iuo or  uo  - s< 


Thus  for  each  s,  W(k,s)  is  cubic  in  k,  with  a  single  root  at  k  =  s 
and  a  double  root  at  k  =  uQ  -  s  =  v.  Since  these  are  the  only  zeros  of 
W,  and  the  leading  coefficient  of  W  is  -2e  <  0,  and  v  <  s,  we  conclude 
that  W  ^  0  for  v  <_  k  <_  s. 

Note  that  in  the  above  proof  e  may  be  arbitrarily  large. 
Furthermore  lira  un  =  e  =  .61  <  1.  Also,  the  double  root  of  W  at 
k  =  v  shows  that  Q  is  a  two  dimensional  analog  of  a  one  dimensional 
scalar  contact  discontinuity. 

Finally,  to  prove  that  T[l,  0:  Uq]  satisfies  the  entropy 
condition  for  all  e  >  0,  it  suffices  to  show 


W(k,s)  -  (2y(s)  -  2sy' (s)) (k-0) 

+  (k2-0)y ' (s)  -  2(k2  +  ek3  -  0)  >  0, 


(4.24) 


for  0  £  k  _<  s,  uQ  =  e  2e  £s  <_  1.  Substituting  (4.3)  and  (4.4)  , 
we  have 

W(k,s)  -  (2(2s  -  2es2ln(s)  +  (e-l)s2)  (4.25) 

-2s  (2  -  4csln(s)  -  2s))  k 
+  k2  (2  -  4£sln(s)  -  2s)  -  2  (k2  +  ek3) 

-  2k  (s2  (2elns  +  e  +  1))  +  k2  (-4esln(s)  -  2s)  -  2ek3. 

Thus  W(k,s)  *  0  if  k  *  0,  or  if  k  =>  s.  We  have  already  seen,  in  (4.5), 

that  ~~  (0,  s)  >  0  if  s  >  Uq,  and  ~  (0,  Uq)  ■  0.  Since  W  is  cubic 
in  k  with  leading  coefficient  -2e  <  0,  we  may  conclude  that  W  >_  0 
for  0  ^  k  £  s,  Uq  <  s  <  1. 
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